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ABSTRACT 

Aims. We present a solution to improve the performance of coronagraphs for the detection of exo-planets. 

Methods. We simulate numerically several kinds of coronagraphic systems, with the aim of evaluating the gain obtained with an 
adaptive hologram. 

Results. The detection limit in flux ratio between a star and a planet (F s /F,,) observed with an apodized Lyot coronagraph charac- 
terized by wavefront bumpiness imperfections of A/20 (resp. /1/100) turns out to be increased by a factor of 10 34 (resp. 10 51 ) when 
equipped with a hologram. 

Conclusions. This technique could provide direct imaging of an exo-Earth at a distance of 1 1 parsec with a 6.5m space telescope such 
as the JWST with the optical quality of the HST 
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r*** 1 . Introduction phase masks could be mitigated w ith devices such as the achro- 

, matic annular groove phase mask (Mawet et al. 20O7h. 

, Most of the 30 or so exo-planets discovered since 1995 , . ... , 

O ■ (lMavoretal.i ri995) have been detected by the radial velocity or ° ne can also ™P rove the exlst ^ n f corona g/ a r* s and era f the 

the photometric transit methods. In a few cases, transits could of , * e stellar point-spread function (PS F) by apodizing 

be observed spectroscopic^ an d provide d data on chemical the pupil, for exampl e using prolate functions j Aime et al. || 2002t 

composition and temperature (Swain et al. 2009). A few sum- I Soummer etalj2003|) . Unfortunately, a large fraction of the light 

H 1 ciently bright planets were imaged fr om the ground, with adap- 18 absorbed by the apodizaton mask. The loss is av oided with an 

! live optics in the ne ar-infrared, fir^v lChauvin et al. ( 2004), and apodizing device ( |Guyon| | 2003t | Guyon et al. || 2007|) using a pair 

more recently by Mar oiTitd] H) . who removed most of the of d f°rted ™™ TS ,° m °t y g dlstnbutlon ac , ross the 

unwanted diffracted starlight by "angular differential imaging". P u P lL W ^ ^ loss " less achromatic apodization exo-planets as 

Another such detection dKalas et al.l& was made, in yellow faint as 10 can ln P rincl P le be dlrectlv lma S ed - 

and red light, with the Hubble Space Telescope. The direct imag- However, in practice, even a theoretically perfect coronagraph is 

ing of weaker extra-solar planets, including exo-Earths about greatly affected by the residual bumpiness of the incoming stel- 

10 10 times fainter than their parent star at visible wavelengths, lar wavefront, caused by imperfect mirror polishing or residual 

remains an enormous challenge since the residual starlight in the atmospheric turbulence, which cannot be perfectly corrected by 

image must be removed very efficiently to allow the detection of adaptive optics. The residual halo of starlight in the image typ- 

faint planet images. ically exhibits a speckle pattern at a relative level much higher 

The coronagra phic te chniques developed since the solar corona- tha n 10 10 Marechal's formula dLabeyrie et al.|2006|) shows that 

graph of iLyoil (fT939l) . particularly in recent years, have signif- t0 detect an exo-planet 10 9 times weaker than its parent star with 

icantly improved its nulling gain, which was limited to about a 10m mirror operating at 1 |^m, the RMS wavefront error among 

10000. The "Lyot coronagraph" w as improved by replacing 10c m patches, such as obtained by 10000 actuators, should re- 

the opaque mask with phase masks dRoddier & Roddieij|l997t ma in below 0.5nm. 

Riaud et al. 20 0ll), whic h attenuated the background level below Since such accuracy is not achievable with the present figuring 
10 -7 dRiaud et al.l l2.Q03). The stronger chromatic dependance of techniques, one method to remove the residues is to add stages 
of adaptiv e elements. Several adaptive devices proposed in the 
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Boursier FRIA 



or adaptive elements, several adaptive devices proposed in me 
literature dCodona & Angelll2004lLabevriell2004l: iPutnam et all 
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Fig. 1. Holographic coronagraph using a modified Lyot train. 
The focal occultor is replaced by a deviating micro-prism which 
separates the star's geometric beam for use as the holographic 
reference beam, while preserving the outer diffracted rings and 
speckles, together with the images of planets, forming the di- 
rect beam. The micro-prism size, amounting in our simulations 
to (8A//Z)) or 13.5 (a.m in visible light at F/3, can vary between 
that of the Airy peak and a dozen rings. Equivalently, the micro- 
prism can also be replaced by a reflective focal plane with a 
hole selecting the reference beam. A larger deviating prism de- 
flects the reference beam in such a way that it intersects the di- 
rect beam in the relayed pupil image where the Lyot stop is lo- 
cated. A dynamic hologram, located at the same position, first 
records the interference fringes produced by both intersecting 
beams, and then reconstructs a copy of the star's direct wave- 
front, without the planet's contribution. The hologram is made 
to introduce a n phase shift in the reconstructed beam, causing 
it to interfere destructively with the "live" stellar wavefront (di- 
rect beam) also transmitted through it by zero-order diffraction, 
thus nulling the star's speckles and improving the planet's con- 
trast in the relayed sky image detected by a CCD camera. The 
lens relays the telescope aperture at the pupil plane and its focal 
plane at the camera. In the final image, the reference beam is 
diffracted by the hologram into several orders (-lref, Oref, and 
+ lref in the figure). The order "+lref" is a reconstruction of the 
speckled wavefront, phase-shifted by n. 



2006) use the following interesting property: if the star is unre- 
solved, the light of the speckle halo is coherent with that of the 
central source, usually absorbed by the Lyot mask. This light in 
principle can be made to interfere destructively with the speckle 
halo. It improves the detectability of the faint non-coherent exo- 
planet, thus relaxing the accuracy of the AO actuators. 

Here, we prese nt simulations w ith the adaptive hologram 
method proposed bv lLabevriel d2004l) . Located close to the Lyot 
stop, a hologram removes most of the residual star light by 
adding to it a phase-shifted copy of its wavefront. In Sect. [2] we 
describe the optical design of the holographic technique adapted 
to the traditional Lyot coronagraphs. In Sect. [3] we discuss the 
adaptive hologram technology. In Sect. [4] we give an analytic 
description of the diffracted orders generated by the hologram 
and discuss how to optimize the parameters of the system. In 
Sect. [5] we show the results of our numerical simulations, dis- 
cussing the gain obtained under different conditions: ideal ones 
(Sect. 15. IT ) and with static wavefront bumpiness and photon noise 
(Sect. |3T2l 15.31 ). Finally, Sect. [6] deals with the degrading effect 
of a star which is partially resolved, its possible mitigation, and 
the design of an achromatized version. 
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Fig. 2. Simulation of an imaged star-planet system in a Lyot 
coronagraph equipped with an adaptive hologram (see Fig. [TJ. 
Here, for clarity, the planet is particularly bright and far from the 
parent star. The diffracted orders -1,0 and + 1 of the planet are 
noted respectively -lpl, Opl and +lpl. The orders -1,0 and +1 
of the star are noted respectively -lstar, Ostar and +lstar. The 
order + lstar looks faint in this figure because it is fainter than 
the direct flux of the star ( Ostar) and the reference beam (Oref). 
The orders - 1 (twin wave), and + 1 of the reference beam are 
noted respectively - lref, Oref and +lref. Star nulling is caused 
by the destructive interference of its order (Ostar) and the or- 
der + 1 of the reference beam (+lref), both superposed on the 
detector. Here, no n phase shift is introduced in the reference 
beam, so as to show the diffracted orders +lref, and Ostar. The 
"twin wave" falls out of the detector at left , but here spuriously 
appears at right due to aliasing, itself caused by the moderate 
sampling which we had to use in the hologram. 



2. Coronagraph design improved with an adaptive 
hologram 

As previousl y described ( Labevrie & Le Corollei] l2004t 
lLabevriel2 004). and shown in Fig.Q] the focal occultor in a Lyot 
coronagraph can be built in the form of a micro-prism deflector 
preserving most star-light to produce an off-set reference wave. 
It interferes with the direct wave in order to generate a hologram 
which is initially recorded in the pupil plane near the Lyot stop, 
and then exploited for nulling most starlight while preserving 
the planet's image formation. Depending on the focal ratios, 
the reference beam may intersect the direct beam at an angle 
of several degrees. We note fa and I d = \ifr d \ 2 the complex 
amplitude of the direct beam and its intensity at a point P of the 
hologram's plane (see Fig. Q]), while ifr r and I r = |^ r | 2 are the 
corresponding values for the reference beam. 
The interference between the reference beam and the direct 
beam creates fringes in the speckles of the hologram (see Figs. [3] 
and|4|. The fringes of the recorded hologram behave like a grat- 
ing, diffracting several orders which become focused in the focal 
plane (see Fig.[2j». We give a more complete and detailed analysis 
of the creation and role of the diffraction orders below. The com- 
plex amplitude resulting from the interference of the two beams 
at point P in the pupil is 



lff=)ff d + xl/ r , 

while the intensity / is 

= I r + ifr d iff* + i/r* d i/f r + I d 



(1) 



I, 



Ir Ir Ir 



(2) 



The hologram has an amplitude transmittance r = F^ 2 whe re y is 
the classical intensity contrast in photographic materials dPerezI 
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Fig. 3. Small part of a hologram obtained with a laboratory sim- 
ulator designed according to Fig[TJ suitably magnified to display 
the low-contrast fringes within the speckles. 



2000). By illuminating the recorded hologram with the reference 
beam only, a reconstructed image of the star's speckles appears 
on the detector, noted +lref on Fig. [2] Using a reference beam 
phase-shifted by n, we obtain the same image phase-shifted by 
7r if the hologram is a positive print. Note that in the article of 
lLabevrie & Le Corollerl (|2004|) the n shift was obtained with a 
negative hologram, which is equivalent. Finally, the order of 
the direct beam (Ostar) adds destructively with the order "+lref", 
thus nulling the residual speckles of the star. 
The planet's light, being incoherent with the reference beam and 
little affected by the micro-prism, is not reconstructed by the 
hologram, and therefore escapes nulling. If we illuminate the 
recorded hologram with the reference beam, now phase-shifted 
by 7T, the transmitted complex amplitude becomes: 



l{l = T (l/v + ll/ r e in ) = T(fff d - lf/ r ) . 



(3) 



If y = 2, the usual value considered optimal in holographic prac- 
tice, the product of both factors having respectively four and two 
terms gives eight terms. As expected, since the reconstructed di- 
rect wave is 7r-shifted with respect to the transmitted direct wave, 
the corresponding pair of terms cancels. Another pair of terms 
also cancels for a related reason. The expression thus simplified 
is: 



i/r = -I ri f, r + I d iff d - itfif/* d + tfr*if/ 2 d . 



(4) 



The remaining terms of Eq.|4]represent various orders diffracted 
by the hologram, locally behaving like a diffraction grating, and 
becoming separated on the camera. Some of them widely spread 
their speckled light on the camera (see Fig. thereby degrad- 
ing the visibility of planet peaks, which are imaged through the 
hologram without being much affected by it. It is of interest to 
discuss the impact of each term on the planet detection: 

1 . the term I r if/ r is the transmitted reference wave, producing an 
intense but sharp focus on the camera (order of the refer- 
ence beam; Oref in Fig . [2j> . Apodizing the reference beam im- 
proves the nulling depth since its Airy rings in its zero-order 
image pattern become attenuated, and thus contaminate less 
the planet's image. Such apodization occurs naturally if the 
micro-prism separator selects the central Airy peak down to 
its first dark ring. Additional spatial filtering can be achieved 
at the micro-prism if needed. 
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Fig. 4. Numerical simulation of the central part of a hologram 
recorded in the pupil with reference beam at 65° incidence. For 
a contrasted display here, the wave bumpiness is adjusted at the 
A/40 level, at A = 550nm, and the reference beam intensity is 
much reduced, to show both the speckles and the finer fringes 
within them. There are approximately four fringes per speckle, 
and four pixels per fringe period. The total pixel count is 1599 x 
1599, only 200 x 200 of which are displayed here. 



the term I^d describes a wave that propagates close to the 
direct stellar residue wave (I^d), nulled by the hologram, 
but with attenuated and spatially modulated intensity. This 
term is called the "weaker wave" in Fig. [2] The modulation 
causes some diffractive spreading on the camera, contami- 
nating the planet image. This is mitigated by using a very 
intense reference beam, relative to the direct beam, so that 
Idftd htyd- The intensity thus achievable, as limited by 
the energy content of the Airy peak in the telescope's focal 
plane, defines the maximal possible gain with a hologram. 
This suggests that the maximum nulling achievable is: 



y i 2 



1 ' 



(5) 



P d 



where the J] p summing extends to all points of the pupil 
plane that are not blocked by the Lyot stop. The numerical 
simulations support this estimate. 

3. The term i/^t^ is the order +1 of the direct beam (+lstar in 
Fig. 0; 

4. ifyf* d is a wave similar to the direct beam (direct stellar 
residue wave), since its complex amplitude is proportional 
to the conjugate of if/ d . This last wave is termed the "twin 
wave" (order -1 of the reference beam; -lref in Fig.|2]l. Its 
contamination is discussed in Sect. [4] 

The order -1 of the direct beam (-lstar) cancels with the wave 
Idftr that propagates close to the order of the reference beam. 

3. Adaptive hologram technology 

A dynamic, or adaptive, hologram works both as a wavefront 
sensor and an actuator array, thus behaving like the feedback 
loop of conventional adaptive optics. In the absence of rewritable 
holographic materials having enough light sensitivity and re- 
sponse speed, the hologram's sensing and playing functions can 
also be achieved by two separate components: a camera and 
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Fig. 5. Approximate S/N ratio of planet detection vs. the 
star/planet flux ratio, in the presence of A/20 wavefront bumpi- 
ness at /I = 550nm and without photon noise, (see Sect. l5.2l ). The 
detection limit (dashed) was chosen at S/N — 3. We find a value 
of Yjtf/Tjld ~ 10 3 for the solutions providing the hologram. 
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Fig. 6. Approximate S/N ratio of planet detection vs. the 
star/planet flux ratio, in the presence of Aj 100 wavefront bumpi- 
ness, at A — 550nm and without photon noise (see Sect. 15. 2b . The 
detection limit (dashed) was chosen at S/N = 3. We find a value 
of Yj ij/T* 1% ~ 10 6 for the solutions providing the hologram. 



a Spatial Light Modulator (SLM). To detect and process the 
light beam simultaneously, the camera can be fed by a beam 
splitter while the SLM is located in the relayed pupil, i.e. at the 
hologram position indicated in Fig. Q] The latter can be driven 
by a video signal from the former, at the image rates of standard 
television, in the absence faster versions. Both components 
are small and commercially available devices which can fit, 
together with the beam-splitter, within a cubic inch. Of interest 
are EM-CCD camera chips, incorporating electron multipli- 
cation which makes them nearly photon-limited at low light 
levels; and SLMs such as those using a liquid crystal film on a 
silicon matrix, currently available with 1280 x 720 pixels (see 
www . cambridgecorrelators . com/products . html). Such 
existing components are readily usable to test a holographic 
coronagraph on a laboratory bench and then on a telescope, 
particularly at red and near infra-red wavelengths where the 
lifetime of speckles is longer. Inserting a computer or dedicated 
fast processor in the video connection is useful to adjust the n 
phase shift, the gamma contrast, etc. . . of the printed hologram. 
Much slower speeds suffice in space, where an artificial star 
such as a remote laser source can in principle be used to record 
the hologram. 

With these existing components, a typical observing sequence 
involves the near simultaneous printing of the hologram with the 
live recorded pattern. 

4. Numerical simulations with optimized hologram 
parameters 

We adjusted the size of the micro-prism to that of the fourth dark 
ring in the Airy pattern; the planet was placed in the fifth ring. 
This reproduces the angular separation between a star and an 
Earth-like planet located at 1 UA from its parent star, at 1 lpc 
from us (like the Sun and Earth as seen from Gliese 436), and 
observed with a 6.5m telescope. In order to simplify the com- 
putations, we coarsely simulated the mirror bumpiness or turbu- 
lence with a square grid of bumps with random amplitudes, and 
verified that the shape of the phase cells on the mirror do not 
affect our results. 



For minimal calculation noise, which represents possible prob- 
lem given the high dynamic range considered, while keeping a 
reasonable computation time and usage of our 2GB RAM mem- 
ory, we used 1599 x 1599 pixel arrays, with a pupil spanning 402 
sampling points. 

The bumpiness of the incoming wavefront is converted into in- 
tensity speckles having random phases in the relayed pupil, i.e. 
in the hologram, in response to the occultation of the Airy peak 
in the focal plane. In the pupil plane, the interference with the 
oblique reference beam generates the finer fringes seen within 
each speckle (Figs. [3] and 0). For a good sampling in the holo- 
gram, four pixels at least are needed per fringe period, and simul- 
taneously at least four fringes per speckle. Increasing the angle 
9 between the direct and reference beams increases the num- 
ber of fringes per speckle, but decreases the number of pixels 
per fringe. It also increases the spacing of the diffracted orders 
in the focal plane, which improves the contrast of the planet's 
peak. In our simulation, we used about 4 pixels per fringe for an 
angle 8 = 65° (see Figs. |4]>. This corresponds to at least 4 fringes 
per speckle, if the typical scale of the bumpiness at the entrance 
pupil is larger than 12 x 12 pixels. For optimal performance, it is 
important to adjust the following parameters: 

- In accordance with the classical theory of holograms, the ref- 
erence beam must be substantially more intense for a faith- 
ful wavefront reconstruction by the hologram, as apparent in 
the second term of Eq. [4] In practice here, the hologram's 
nulling gain is proportional to the square intensity ratio, i.e. 
as 10 3 in our numerical simulation with wavefront bumpi- 
ness imperfections at A/20 at A = 550nm. This ratio is lim- 
ited in practice by the energy content of the Airy peak, di- 
verted by the micro-prism to form the reference beam. The 
size of the micro-prism or attached filtering aperture, possi- 
bly apodized, should be optimized for maximal light collec- 
tion while keeping a flat wave front. We found that an aper- 
ture sized like the Airy peak is optimal in this respect and 
also to minimize the inner working angle. 

- As in conventional Lyot coronagraphs, the size of the Lyot 
stop was adjusted to block the bright edge of the pupil image. 
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Table 1. Summary of the approximate limiting detection of the 
flux ratio F s /F p with the different configurations mentioned in 
Sect. [5] without photon noise 1 . 



- Based upon Eq. [4] we conclude that an optimal value for y 
is 2 using a positive print, with a phase shifter then intro- 
duced in the reference beam during read-out, and -2 if it is 
negative. 

- The twin wave is described by the term \jA in Eq.|4] It gives 
a focal image identical to the direct wave, but shifted, and 
non-overlapping if the reference beam angle, with respect to 
the direct beam, exceeds the apparent pupil size. In our sim- 
ulations, however, the pixel sampling of the hologram was 
insufficient to properly generate the twin wave, causing it to 
be aliased by the Fourier transform algorithm, and to appear 
in the final image at an incorrect "folded" location which 
contaminated the planet image (See Fig.|2]i. The aliasing ef- 
fect causes the twin wave to be nearer to the order of the 
direct beam (and then to disturb the flux of the planet) than 
if it was not aliased. 

In order to remove such numerical effects, we simulated a 
suppressed twin wave by removing the corresponding term, 
in complex amplitude, in the numerical calculation, and 
found an improved detection sensitivity, part of which may 
result from the removed aliasing effect. 
In this paper, we did not study the physical effect of the twin 
wave, but we know that the solution is intermediate between 
the case including the twin wave (C in Table. [TJ and the case 
where its effect is analytically subtracted (D in Table. [TJ. 
More realistic simulations of the twin wave's effect are de- 
sirable in furthe r work, and this may r equire non-FFT calcu- 
lation methods ( Soummer et al . 2003). 
Nevertheless, twin waves can be eliminated by using a thick 
hologram, also known as a Lippmann-Bragg hologram, in 
order to work in the case (D) of Table. [JJ These have strat- 
ified nodal planes, rather than fringes, which selectively re- 
flect a single first-order wave. Then, the hologram works like 
a "blazed grating", sending most light in the reconstructed 
image (order +lref). 



5. Gain evaluation 

Using the optimizations listed in Sect. [4] we were able to re- 
produce a wide range of coronagraph configurations, in order to 
compare their planet detection limit versus the "star/planet" flux 
ratio (F s /F p ). 

(A) Classical Lyot coronagraph; 

(B) Apodized Lyot coronagraph; 

(C) Apodized Lyot coronagraph with hologram and an apodized 
reference beam; 

(D) Apodized Lyot coronagraph with hologram, apodized refer- 
ence beam, and subtracted twin wave. 
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Fig. 7. The four rows show simulated images with the A, B, C 
and D coronagraph types mentioned in Sect. [5] The star / planet 
flux ratio is 10 3 in the left column, 10 5 5 in the central column, 
and 10 6 5 in the right column. The wavefront bumpiness is A/20 
as described in Sect. 15.21 The intensity scales at right indicate 
the number of photons per pixel, for a star with magnitude my = 
7, a telescope diameter D = 6.5m, A = 550nm and spectral 
bandwidth AA = lOnm. The exposure time is 60s for both the 
hologram and the science camera 



The performance is evaluated under three different conditions: 
ideal ones (no mirror imperfections and absence of photon 
noise); with the introduction of mirror imperfections and finally 
with photon noise in addition to the mirror imperfection. These 
steps are treated in the following subsections, and summarized 
in Table [TJ for the two first conditions. The speckle noise was 
coarsely evaluated in the images using the following equation: 



I peak I speckle 



speckle 



(6) 



1 "/1/100" and "/i/20" are referred to the gain after the introduction 
of wavefront bumpiness imperfection at A/ 100 and A/20 (Sect. [5"2t . 



where I pea k is the mean intensity at the position of the planet; 
I speckle is the mean intensity of the speckles near the planet posi- 
tion, and cr speckle is the corresponding root mean square fluctua- 
tion. 



5.1. Gain with a perfect mirror 

The first results are obtained under the assumption of perfect 
conditions, in the absence of mirror bumpiness and photon noise. 
For the classical Lyot coronagraph (A), we obtain a detec- 
tion limit in flux ratio of F s /F p = 10 3 2 . In the case of the 
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Star flux / Planet flux 

Fig. 8. Approximate S/N ratio as a function of the star/planet 
flux ratio, evaluated with /i/20 wavefront bumpiness and photon 
noise (see Sect. 15. 3t for the coronagraph (D). curves are calcu- 
lated for different stellar magnitudes, assuming D = 6.5m, 60s 
exposures, A = 550nm and AA = lOnm. The performances of the 
apodized Lyot coronagraph at my — 7 are also plotted (squares). 
The detection limit (dashed line) was chosen at S/N = 3. 



apodized Lyot coronagraph (B), the detection limit is increased 
to 10 110 . Introducing the hologram in the optical scheme with 
the apodization of the reference beam (C), the limit is 10 11 °. 
Our simulations appear to be limited to this range by numerical 
noise. The hologram becomes most valuable in the presence of 
mirror imperfections. 

5.2. Gain with bumpy mirror 

In order to test the performance under more realistic conditions, 
we have introduced a random static bumpiness on the incoming 
wavefront, with a A/20 peak-to-valley amplitude. We also per- 
formed simulations at Aj 100, a situation intermediate between 
perfect conditions and A/ 20 amplitude. The results are shown in 
Figs.HtoH] 

The detection limit in flux ratio with a A/ 20 wavefront bumpi- 
ness is F s /F p - 10 3 5 with the classical Lyot coronagraph (A); 
10 4 - 5 with the apodized Lyot coronagraph (B); 10 6 2 with the 
hologram and apodized reference beam (C); and 10 7 9 after sub- 
tracting the twin wave (D). At A/100 we find FJF p = 10 3 5 with 
the classical Lyot (A); 10 5 7 with the apodized solution (B); 10 76 
with the hologram and apodized reference beam (C); and 10 ll) 8 
after subtracting the twin wave (D). 

The results with mirror bumpiness reveal the large gain then ob- 
tained with the introduction of the hologram, especially when 
subtracting the twin wave. Indeed, with A/20 wavefront bumpi- 
ness, it is not possible to detect a planet with an apodized Lyot 
coronagraph if the ratio F s /F p is larger than 10 4 5 (see Fig. [7j>. 
This limit is pushed to 10 79 using the same coronagraph 
equipped with a hologram and after the subtraction of the twin 
wave. The detection limit is increased by a factor 10 3 4 . This 
number is about equal to the intensity ratio, predicted by the an- 
alytical estimation (see Eq.|5]). 



5.3. Gain with bumpy mirror and photon noise 

After having tested the coronagraphs under perfect conditions 
and in the presence of mirror bumpiness, we added photon noise. 
Unless the hologram is recorded with infinitely many photons, 
its recorded fringes are noisy, and this degrades the nulling depth 
in the camera image by creating a broad speckled halo, which 
degrades the visibility of the planet's peak. It can be shown that 
the sensitivity limitation, regarding the detection of faint plan- 
ets, is then ultimately the same as with a conventional adap- 
tive system feeding a perfect coronagraph or apodizer, where 
the nulling depth is similarly limited by the number of pho- 
tons detected by the wave se nsor. The calculation given in the 
former case dLabevriel 12004) indeed also applies to the latter: 
N p photons detected by either a hologram or a conventional 
wave sensor limit the achievable peak/halo nulling depth G to 
Gmax = N p . In b oth cases, as well as with a Mach-Zehnder in- 
terference nuller dCodona & Angefl2004t|Putnam et al.ll2.Q06l) . it 
follows that recording and observing exposures lasting the same 
time leave at least one star photon per speckle in the nulled cam- 
era image, which can be low enough to detect planet peaks con- 
taining several photons. The holographic nuller has no theoreti- 
cal advantage in this respect. 

Hybrid forms using both conventional adaptive optics for coarse 
wavefront correction, and a hologram for fine correction are also 
possible. They can use a single wave sensor such as a camera 
located in the hologram plane and serving as the holographic de- 
tector. Its signal can indeed provide a wavefront map to activate 
a deformable mirror in the entrance aperture. It can also be fed to 
a dynamic holographic plate, incorporating its own actuators, in 
the form of the fringe patterns acting as tiny gratings within each 
of its speckles. Like ordinary actuators in a servo loop, their per- 
formance is degraded by photon noise in the control signals, i.e. 
the recorded fringes with their "frozen" photon noise. During the 
observation, after the hologram recording stage, the more usual 
form of "live" photon noise is also present in the recorded im- 
age and further contributes to degrading the detectability of a 
planet's peak if its level is not much higher than the surrounding 
speckle peaks. 

In our simulations, we have included both the "frozen" and the 
"live" photon noise contributions. We assumed that exposures of 
equal duration served for the two phases of hologram recording 
and observation. This may be optimal if the hologram is recorded 
with the same star which is subsequently observed. If, however, 
a brighter reference star or a laser artificial star serves to record 
the hologram, the results given below may be considered as a 
lower limit since the photon noise of the hologram becomes at- 
tenuated in this latter case. Poisson-distributed photon noise was 
generated in the hologram an d the final image u sing a standard 
IDL routine, based on C code dPress et al.lll992l) . 

Starting from the best solution, i.e. a hologram with an apodized 
reference beam and without a twin wave, we assumed 60s ex- 
posures, both in recording the hologram and observing, using 
a 6.5m aperture in the V band (550nm), with a filter width of 
lOnm. The results are shown in Figs. [8] and [9] (simulations at 
A/20 and Aj 100 respectively). 

We see that the photon noise is negligible for stellar magnitudes 
brighter than 7 (at A/20) and 3 (at Aj 100); the hologram perfor- 
mance is then maximal, and independent from the star magni- 
tude. In our simulation, those magnitudes correspond to about 
« 3 x 10 6 (at ,1/20) and * 3 x 10 s (at A/100) photons per speckle 
in the hologram, and their performances are degraded below this 
value. 
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Performances while adding photon noise 

Hoi, apodized ref, without twin wave 
Wavefront bumpiness: A./100 




10° 10 2 10 + 10 5 10 s 10 la 10 12 
Star flux / Planet flux 

Fig. 9. Approximate S/N ratio as a function of the star/planet 
flux ratio, evaluated with Aj 100 wavefront bumpiness and pho- 
ton noise (see Sect. 15. 31 ) for the coronagraph (D). curves are cal- 
culated for different stellar magnitudes, assuming D = 6.5m, 60s 
exposures, A = 550nm and AA = lOnm. The performances of the 
apodized Lyot coronagraph at my — 7 are also plotted (squares). 
The detection limit (dashed line) was chosen at S/N = 3. 

Above magnitude 13 (at A/20) and 11 (at X/ 100), where the 
hologram has ~ 3 x 10 4 (at A/20) and « 2 x 10 5 (at A/100) 
photons per speckle, the sensitivity gain vanishes with respect to 
a Lyot apodized coronagraph without a hologram. 
With a ground-based telescope, affected by "seeing", bright stars 
may provide enough photons during brief exposures, shorter 
than the lifetime of "seeing", to activate both a conventional 
adaptive optics system giving a high Strehl ratio, and, within 
the coronagraphic attachment, the recording of a dynamic holo- 
gram. On fainter stars, a bright Laser Guide Star may similarly 
serve for both stages of adaptive correction. 
In space, where the wavefront bumpiness is greatly reduced and 
varies much more slowly, some planet detection sensitivity can 
be gained by recording the hologram on a star brighter than the 
star observed, or on a ground-based laser source (which is not 
affected by turbulence if its emitting aperture is unresolved). 

It could be of interest if the hologram could attenuate the 
fixed stellar speckles escaping the adaptive correction. We at- 
tempted to address this by replacing ^ with a sum of a constant 
and a randomly variable term, in Eq. [2] but did not elucidate the 
matter. Further work with simulations will be useful. 

6. Discussion 

In the following subsections, two aspects for the future devel- 
opment of the instrument are considered: the chromatism of the 
coronagraph with a hologram; and the problems introduced by a 
star that is poorly resolved. 

6.1. Achromatizing a hologram 

The fringe spacing in the hologram is normally proportional to 
wavelength but can be made invariant if the deviating prism 
(See Fig. [TJ in the reference beam is replaced by a diffrac- 
tion grating. If the grating operates in the first order, its angu- 
lar dispersion indeed increases the incidence angle of light at 
increasing wavelengths. This can reduce the hologram's wave- 



length sensitivity, both during the recording stage and the ob- 
servations, and increase the usable spectral bandwidth, although 
the speckles contributed by the direct beam are also wavelength- 
dependant. Further simulations, achieved with a range of wave- 
lengths, would be of interest to specify the bandwidth then 
achievable. 

Another approach involves a Lippmann-Bragg hologram, 
i.e. a thick hologram where the fringes are patterned as strat- 
ifications throughout the depth of the recording layer. These 
can be wavelength-multiplexed, and can simultaneously recon- 
struct the variously colored wavefronts that have been recorded. 
Some recording materials, such as lithium niobate single crystals 
or polymeric compounds (Shishido & Ikeda 2007), are erasable 
and re-usable. Their moderate recording sensitivity however 
may require a bright laser star. 

6.2. Resolved parent star 

The original Lyot coronagraph is highly tolerant of a star be- 
ing resolved, as demonstrated by its initial success on the solar 
corona. However, a star which is slightly resolved by the tele- 
scope or hypertelescope dLabevriel 1996t lLe Coroller et al.l2 004) 
can degrade the hologram recorded with its light, and also de- 
grade the nulling depth in the image cleaned by a pre-recorded 
hologram. In both cases, what matters is the invariance of the 
hologram, particularly in terms of fringe positions, with respect 
to a slight motion of a point star. Such motion produces iden- 
tical translation shifts of the central and peripheral focal pat- 
terns separated by the micro-prism. The fringe pattern recorded 
in the hologram, located at a nearly infinite distance in the re- 
layed pupil, is therefore nearly invariant. The pupil indeed re- 
mains fixed, and both interfering wavefronts reaching it become 
tilted by the same amount, while their phase difference is in- 
variant. Their speckle detail is however slightly modified by the 
varying edge effects at the boundaries of the micro-prism and 
field lens. The hologram is thus expected to be somewhat toler- 
ant of a resolved star serving to record it. Similarly, a hologram 
recorded on a point source and then used to observe a resolved 
star is also expected to efficiently null its coronagraphic residue, 
since the live fringes remain contrasted. Further simulations and 
laboratory experiments will be of interest to estimate the tolera- 
ble apparent size of the star vs. the desired nulling depth. 

7. Conclusions 

Our diffractive analysis and simulations of a holographic coro- 
nagraph, and the comparison with other forms of adaptive coro- 
nagraphy, show that the theoretical photon-limited sensitivity in 
detecting faint exo-planets is comparable. These methods will 
therefore need to be compared in terms of their technical imple- 
mentations. Hybrid methods, for example combining an adap- 
tive mirror before the occultor and a hologram after, are likely 
to be similarly limited by the photon noise but may also have 
practical merits, suc h as relaxing the accuracy of the actuators 
dPutnam et ai1l2006l) . 

Depending on the holographic processes and materials that will 
become available, various forms of practical implementation 
may be of interest. With some re-writable holograms using am- 
plifying photo-sensitive materials, no electronic image process- 
ing may be needed. Instead, and for a better flexibility, the 
recording and diffractive functions may be separated, using re- 
spectively a camera and a spatial light modulator. 

With A/ 20 wavefront bumpiness, the planet detection limit is 
improved, in terms of star/planet flux ratio, from 10 4 5 to 10 7 9 . In 
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terms of photon noise, we found that the hologram improves the 
performance of the coronagraph if it is recorded with more than 
~ 3 x 10 4 photons per speckle. With A/ 100 wave front bumpi- 
ness, the detection limit improves from 10" to 10 10 - 8 , and the 
hologram improves the performance of the coronagraph if it is 
recorded with more than « 2 x 10 5 photons per speckle. 
We also discussed ways of making this holographic technique 
achromatic, and its tolerance for a poorly resolved star. In or- 
der to reach the best results, a thick Lippmann-Bragg hologram 
is of interest to remove the twin wave, but it restricts the range 
of available holographic components. Following laboratory sim- 
ulations, already initiated, tests on ground-based adaptive tele- 
scopes may be possible with available camera and Spatial Light 
Modulator components. 
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